Exploring the Universe beyond the Photon Window* 

Luis A. Anchordoqui 

Department of Physics, Northeastern University, Boston MA 02115 



In this talk I review how to identify cosmic ray accelerators that are high 
energy neutrino emitters. I also delineate the prospects for a new multi- 
particle astronomy: neutrons as directional pointers + antineutrinos as 
inheritors of directionality. 

PACS numbers: 96.40.-z, 95.85.Ry, 98.70.Sa 

1. Introduction 

Conventional astronomy spans about 18 decades in photon wavelengths, 
from 10 cm radio- waves to 10 -14 cm 7-rays of GeV energy. Because the uni- 
verse is opaque to photons of TeV energy and above (see Fig. [J) 1 , present 
studies focus on hadrons, neutrinos, and gravitational waves as messengers 
probing the high energy universe. The best candidates to serve as messen- 
gers in a new astronomy of the high energy behavior of distant sources are 
neutral particles. This is because the orbit of a charged cosmic ray can 
be substantially bent, both by extragalactic magnetic fields [2] and by the 
ambient magnetic field of our own Galaxy [3]; destroying the possibility of 
locating the source. The most promising messenger is the neutrino: it can 
be copiously produced in cosmic beam dumps and can traverse unscathed 
dense astrophysical environments. In this talk I delineate the prospects to 
identify high energy neutrino emitters. 

2. Cosmic Ray Astronomy 

For primary energy > 1 GeV, the observed cosmic ray intensity can be 
described by a series of power laws with the flux falling about 3 orders of 
magnitude for each decade increase in energy [I]. In recent years, a some- 
what confused picture vis-a-vis the energy spectrum and arrival direction 
distribution has been emerging. Since 1998, the AGASA Collaboration has 
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Fig. 1. Left: Mean interaction length for photons on the ultraviolet, visible, in- 
frared, and microwave backgrounds £Q. Right: Upper end of the cosmic ray energy 
spectrum as observed by 5 different experiments (for details see main text). 



consistently reported 5 a continuation of the spectrum beyond the ex- 
pected Greisen-Zatsepin-Kuzmin (GZK) cutoff which should arise at 
about 10 10 ' 9 GeV if cosmic ray sources are at cosmological distances. This 
theoretical feature of the spectrum is mainly a consequence of interactions 
of the primary cosmic ray with the microwave background radiation. In 
contrast, the most recent results from HiRes j2J describe a spectrum which 
is consistent with the expected GZK feature. The discrepancy between the 
2 estimated fluxes is shown in Fig. ^ (for comparison, the intensity as seen 
by the Fly's Eye 0, the Haverah Park 9 and the SUGAR [TO] experi- 
ments is also shown). This situation is suggestive of the challenge posed by 
systematic errors in these types of measurements. Further confusing the is- 
sue, the AGASA Collaboration reports observations of event clusters which 
have a chance probability smaller than 1% to arise from a random distribu- 
tion [IJ, whereas the recent analysis reported by the HiRes Collaboration 
showed that their data are consistent with no clustering among the highest 
energy events ^21- I n m Y opinion, it is very important to rigorously define 
the corresponding budget of statistical significance and search criteria prior 
to studying the data, since defining them a posteriori may inadvertently 
introduce an indetermimant number of "trials" and thus make it impossi- 
ble to assign the correct statistical significance to the search result. In this 
direction, with the aim of avoiding accidental bias on the number of trials 
performed in selecting the angular bin, the original claim of the AGASA 
Collaboration was re-examined considering only those events observed after 
the original claim fTHJ. This study showed that the evidence for clustering 
in the AGASA data set is weaker than was previously supposed, and is 
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consistent with the hypothesis of isotropically distributed arrival directions. 
The confusing experimental situation regarding the GZK feature should be 
resolved in the near future by the Pierre Auger Observatory ^1] , which will 
provide not only a data set of unprecedented size, but also the machinery 
for controlling some of the more problematic systematic uncertainties. 

With the controversy over the GZK cutoff, one may have missed the fact 
that below 10 10 GeV AGASA has revealed a correlation of the arrival direc- 
tion of the cosmic rays to the Galactic Plane (GP) at the 4cr level |15j . The 
GP excess, which is roughly 4% of the diffuse flux, is mostly concentrated in 
the direction of the Cygnus region, with a second spot towards the Galactic 
Center (GC). The anisotropy signal spans a narrow energy window, from 
10 8 ' 7 GeV to 10 9 ' 3 GeV. Evidence at the 3.2a level for GP enhancement in 
a similar energy range has also been reported by the Fly's Eye Collabora- 
tion |16j . Interestingly, the full Fly's Eye data include a directional signal 
from the Cygnus region which was somewhat lost in an unsuccessful at- 
tempt to relate it to 7~ray emission from Cygnus X-3 JZ]- 1 Additionally, 
a 3.4(7 excess from an extended region surrounding Cygnus X-3 has been 
reported by the Yakutsk Collaboration JH]- On the other hand, the exis- 
tence of a point-like excess in the direction of the GC has been confirmed 
via independent analysis |20| of data collected with the SUGAR facility. 
This is a remarkable level of agreement among experiment using a variety 
of techniques. 

Independent evidence may be emerging for a cosmic accelerator in the 
Cygnus spiral arm. The HEGRA experiment has detected an extended TeV 
7-ray source in the Cygnus region with no clear counterpart and a spectrum, 



not easily accommodated with synchrotron radiation by electrons [2J- Es- 
pecially intriguing is the possible association of this source with Cygnus- 
OB2, a cluster of more than 2700 (identified) young, hot stars with a total 
mass of ~ 10 4 solar masses [22]. Proton acceleration to explain the TeV pho- 
ton signal requires only 0.1% efficiency for the conversion of the energy in 
the stellar wind into cosmic ray acceleration. Also, the stars in Cygnus-OB2 
could be the origin of time-correlated, clustered supernova remnants form- 
ing a source of cosmic ray nuclei. An immediate consequence of this picture 
is the creation of free neutrons via nuclei photodisintegration on background 
photon fields. These liberated neutrons are presumably responsible for the 
observed directional signals. This implies that it may not be a coincidence 

The complete isotropy at PeV energies revealed by KASCADE data |18f vitiate 
direction-preserving photons as primaries. 
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that the signal appears first at energies where the neutron lifetime allows 
propagation distances of galactic scales, i.e., 10 kpc. 



Neutrinos can serve as unique astronomical messengers. Except for os- 
cillations induced by transit in a vacuum higgs field, neutrinos propagate 
without interactions between source and Earth, providing powerful probes of 
high energy astrophysics (2H1- The deployment of under-ice/water telescopes 
in the Northern |24| and Southern [23] hemispheres will greatly increase the 
statistics required for the realization of such a program. In this section we 
study possible mechanisms for neutrino production in astrophysical sources 
and discuss how flavor oscillations can distort the initial injection spectra. 



In recent years, stronger and stronger experimental evidence for neutrino 
oscillations has been accumulating. Neutrino flavor change implies: (i) that 
neutrinos have nonzero masses, i.e., there is a spectrum of 3 or more neutrino 
mass eigenstates that are analogues of the charged lepton mass eigenstate l a 
(a = e, /i,r); (ii) leptonic mixing, i.e., the weak interaction coupling the W 
boson to a charged lepton and a neutrino, can also couple any charged lepton 
mass eigenstate l a to any neutrino mass eigenstate Uj (J = 1, 2, 3, . . .). The 
superposition of neutrino mass eigenstates produced in association with the 
charged lepton of flavor a, 



is the state we refer to as the neutrino of flavor alpha, where U*j are ele- 
ments of the neutrino mass-to-flavor mixing matrix, fundamental to particle 
physics. 

Throughout this talk we consider 3 neutrino species. In this case, atmo- 
spheric data indicate that i/„ and v T are maximally mixed and reactor 
data [27] points to | C^ e 3 1 2 -C 1. Thus, to simplify the discussion hereafter 
we use the fact that |?7 e 3| 2 is nearly zero to ignore CP violation and assume 
real matrix elements. With this in mind, one can define a mass basis as 
follows, 



3. Neutrino Astronomy 



3.1. Flavor metamorphosis 
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where 9q is the solar mixing angle and 

IO = ^(M-K» (6) 

is the eigenstate orthogonal to \v^). Inversion of the neutrino mass-to-flavor 
mixing matrix leads leads to 

\v e ) = cos8 Q \vi) — sin 00 1 1/2) (7) 

and 

I v*) = sin 6q + cos 0®\v2) ■ (8) 
Finally, by adding Eqs. (JSJ) and © one obtains the flavor eigenstate, 

= —j= [|i/ 3 ) + sin0 |i/i) + cos9 Q \u 2 )] , (9) 

and by substracting these same equations the v T eigenstate. 

The evolution in time of the ut component of a neutrino initially born 
as v a in the rest frame of that component is described by Schrodinger's 
equation, 

|^(r i )) = e-™|^(0)) , (10) 

where rri{ is the mass of i/j and T{ is the proper time. In the lab frame, 
the Lorentz invariant phase factor may be written as e~^ Eit ~ PiL ' , where t, 
L, Ei, and pi, are respectively, the time, the position, the energy, and the 
momentum of z/j in the lab frame. Since the neutrino is extremely relativistic 



t ~ L and Ei = <J p 2 + rn 2 ~ p + mf/2p. Hence, from Eq. © it follows that 
the state vector of a neutrino born as v a after propagation of distance L 
becomes 

\v a {L))^Y, U ^ e ~ l(mV2E)L \^) > ( n ) 

i 

where E k, pis the average energy of the various mass eigenstate components 
of the neutrino. Using the unitarity of U to invert Eq. from Eq. 
one finds that 



/3 



/8» 



^) • (12) 



In other words, the propagating mass eigenstates acquire relative phases 
giving rise to flavor oscillations. Thus, after traveling a distance L an initial 



G 
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state v a becomes a superposition of all flavors, with probability of transition 
to flavor (3, \{vf3\v a (L))\ 2 , given by 

P{y a -> vp) = 5 a/3 -A^2U ai Upi U aj Upj sin 2 Ay , (13) 

i>j 

where Ay ~ 5mf : jL/2E, and 5m 2 j = rrii — rrij. 

For Ay 3> 1, the phases will be erased by uncertainties in L and E. 
Consequently, averaging over sin 2 Ay one finds 

P{v a -» = <5 a/3 - 2 £ [/^ C/ aj [fo . (14) 
Now, using 2£a>j = Sij ~~ Z)t=j> Eq. lfTl|) can be re-written as 

P{v a — > I/g) = — 2J f^m L^/3i ^oj ^/9j + X! ^ ^ ^ 

= <w - + £ ■ ( 15 ) 

Since 5 Q/ 3 = (5 2 ^, the first and second terms in Eq. lfTo|l cancel each other, 
yielding 

P{u a ^u p )=Y J U 2 ai Ul i . (16) 

i 

The probabilities for flavor oscillation are then 

Piy^ = i/ T ) = - (cos 4 # + sin 4 9 Q + 1) , (17) 

■P(f/x ->• f e ) = P{v e -> f M ) = P(^ e -> i/ r ) = sin 2 8q cos 2 6> , (18) 

and 

P(u e -> i/ e ) = cos 4 6> + sin 4 6> . (19) 

Now, let the ratios of neutrino flavors at production in the cosmic sources 
be written as w e : w^ : w T with J2a w a = 1, so that the relative fluxes of 
each mass eigenstate are given by Wj = J^a^a U^j- From our previous 
discussion, we conclude that the probability of measuring on Earth a flavor 
a is given by 

P Va detected = £ £ Wf, U% . (20) 

3 P 

Straightforward calculation shows that any initial flavor ratio that contains 
w e = 1/3 will arrive at Earth with equipartition on the three flavors. 
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3.2. The j-v connection 

There are two principal mechanisms for TeV gamma ray production: 
(i) Electrons undergo bremsstrahlung in the magnetic field and/or inverse 
Compton scattering in the ambient photon sea or (ii) the gamma rays are 
directly traced to 7r° decay. 2 Only the second scenario can accommodate 
baryonic cosmic ray production. Since such cosmic rays are observed, it is 
reasonable to assume that at least some gamma ray sources operate accord- 
ing to the second mechanism. 

Inelastic pp collisions lead to roughly equal numbers of 7r°'s, 7r + 's, and 
7r~'s, hence one expects two photons, two z/ e 's, and four v^'s per ir°. On 
average, the photons carry one-half of the energy of the pion. The average 
neutrino energy from the direct pion decay is (E^) 71 " = (1 — r)E n /2 ~ 
0.22 E n and that of the muon is (E^) n = (1 + r) EJ2 ~ 0.78 E n , where r is 
the ratio of muon to the pion mass squared. Now, taking the from muon 
decay to have 1/3 the energy of the muon, the average energy of the v„ 
from muon decay is {E U/i ) >1 = (1 + r)E n /Q = 0.26 E w . This gives a total 
energy per charged pion (E„ ) ~ 0.48 E n , with a total (£'^ I ) total = 0.96(£' 7 ) 
for each triplet of ir + , ir~ , and 7r° produced. For simplicity, hereafter we 
consider that all neutrinos carry one-quarter of the energy of the pion. 

The total number of 7-rays in the energy interval (E\/2, E%/2) is equal 
to the total number of charged pions in the interval (Ei, E?) and twice the 
number of neutral pions in the same energy interval, 



£1/2 



-dEr, 



dE 7 



-dE n = 2N n o . 



(21) 



Additionally, since N n ± = 
(Ei /A, E 2 /4) scales as 

"£2/4 dF v 



2N w o, the number of in the energy interval 



E-l/4 



dE v 



-dE u 



E * dFl 



dE % 



-dE^ = 2N n 



Now, taking d/dE2 on each side of Eqs. (f2"T|) and l|22|) leads to 



1 dE. 



2 dE~, 



dF2 



--E2/2 



dE w 



and 



\dF v 



E 2 



4 dE v 



dF± 



--E2/4 



dE^ 



(22) 



(23) 



E 2 



respectively. The energy-bins dE scale with these fractions, and we arrive 
at 



dE, 



dF^ 



=E n /2 



dEir 



2 Pions can be produced in pp and/or P7 collisions. In this talk we focus attention on 
the former; a generalization to the photopion production process is straightforward, 



sec 



for details. 
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dF Va 



dE„ 



E v =E n /i 



16 



E-k 



(24) 



dEL 



Err 



for the total fluxes at the source, where ir denotes any one of the three pion 
charge-states. In propagation to Earth a distance longer than all oscillation 
lengths, flavor changing amplitudes are replaced by probabilities. Using 
Eqs. Q17|). (|18|). and ()19|) one can check that the initial flavor ratio 1:2:0 
mutates in a nearly identical flux for each of the three neutrino flavors which 
is equal to [55] 

dF v 



dE v 



2 dF, 



Eu — E'y 1 2 



(25) 



IceCube is, perhaps, the most promising route for neutrino detection 12 
This telescope will consist of 80 kilometer-length strings, each instrumented 
with 60 10- inch photo- multipliers spaced by 1.7 m. The deepest module 
is 2.4 km below the ice surface. The strings are arranged at the apexes 
of equilateral triangles 125 m on a side. The instrumented detector vol- 
ume is a cubic kilometer. A surface air shower detector, IceTop, consisting 
of 160 Auger-style Cerenkov detectors deployed over A e fi ~ 1 km 2 above 
IceCube, augments the deep-ice component by providing a tool for calibra- 
tion, background rejection and air-shower physics. Muons can be observed 
from 10 2 GeV to 10 9 GeV. Cascades, generated by u e , V e , u T , and V T can 
be observed above 10 2 GeV and reconstructed at energies somewhat above 
10 4 GeV. For z^'s of TeV energy, the angular resolution ~ 0.7° allows a 
point source search window of r2i° x i° ~ 3 x 10 -4 sr. 

The Crab Nebula is generally taken as the standard candle of steady 
TeV 7-ray emission. In the energy 1 TeV < £L < 20 TeV, the Crab data 
can be fitted by a single power law [3*0] . 



~dK 



2.79 



±0.02 
±0.5 



x 10" 



Ej 
TeV 



-2.59 



±0.03 
±0.05 



111 



TeV" 



(26) 



Even though the energy spectrum of the Crab Nebula as measured by 
the HEGRA system is in very good agreement with calculations of in- 
verse Compton scattering [SJ, it is interesting and desirable to have an 
independent observational discriminator between these two scenarios. For 



v — 



1 TeV, the number of v., 



showers expected to be detected at 



IceCube is given by 
dN 



A 



elf 



signal 



dF 

dE v -^{E y ) p(E v 



(27) 
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where p(E u ) « 1.3 x 1CT 6 (E u /TeV) - 8 denotes the probability ( generic to 
ice/water detectors) that a v (or v) with energy E u on a trajectory through 
the detector produces a signal [23]. Using Eq. Q25JI straightforward cal- 
culation shows that diV/ eft | signal ~ 11 yr _1 32 . The event rate of the 
atmospheric ^-background is 

— =A cS [ dE v J U+P {E„) p{E v ) A^io xl o « 1.5 yr- 1 , 

«* background JB™" 1 

(28) 

where J u +p(E u ) is the + f M atmospheric flux in the direction of the Crab 
region (about 22° below the horizon) [33]. We conclude that the neutrino 
signal can be easily isolated from background at IceCube. Therefore, the 
hadronic nature of the high energy emission from the Crab Nebula can be 
confirmed or disproved in a few years of operation. 



3.3. The n~v connection 

Cosmic ray experiments have identified an excess from the region of 
the GP in a limited energy range around 10 9 GeV. This is very suggestive 
of neutrons as candidate primaries, because the directional signal requires 
relatively-stable neutral primaries, and time-dilated neutrons can reach the 
Earth from typical Galactic distances when the neutron energy exceeds 
10 9 GeV. In what follows we show that if the Galactic messengers are neu- 
trons, then those with energies below an 10 9 GeV will decay in flight, pro- 
viding a flux of cosmic antineutrinos above a TeV which would be observable 
at IceCube. 

The basic formula that relates the neutron flux at the source (dF n /dE n ) 
to the antineutrino flux observed at Earth (dF^/dE-^) is 



(i_ e -&) /V^m 



dEy J dE n V J Jo de 

1,1 d COS 0J7 



5 

-i 2 



Ep-E n e v {l + cos dv)/m n . (29) 



The variables appearing in Eq. (|29jl are the antineutrino and neutron en- 
ergies in the lab (Ep- and E n ), the antineutrino angle with respect to the 
direction of the neutron momentum, in the neutron rest-frame (#77), and 
the antineutrino energy in the neutron rest-frame (e-jj). The last three 
variables are not observed by a laboratory neutrino-detector, and so are 
integrated over. The observable £77 is held fixed. The delta-function re- 
lates the neutrino energy in the lab to the three integration variables. The 
parameters appearing in Eq. (|29|) are the neutron mass and rest-frame life- 
time (m n and r n ). Finally, dP/de-p is the normalized probability that the 
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decaying neutron produces a V with energy ep- in the neutron rest-frame. 
Note that the maximum V energy in the neutron rest frame is very nearly 
Q = m n — m p — m e = 0.71 MeV. Integration of Eq. (|29j) can be eas- 
ily accomplished, especially when two good approximations are applied. 
The first approximation is to think of the /3-decay as a 1 —* 2 process of 
bvn n — > e~ + V, in which the neutrino is produced monoenergetically in 
the rest frame, with ep = eo — 5m]y(l — m 2 /5 2 miy)/2 ~ 0.55 MeV, where 
<5mjv — 1-30 MeV is the neutron-proton mass difference. In the lab, the 
ratio of the maximum V energy to the neutron energy is 2eo/m n ~ 10 -3 , 
and so the boosted V has a spectrum with Ep- S (0, 10 _3 £^ n ). The second 
approximation is to replace the neutron decay probability 1 — e~ Dm "^ E " Tn 
with a step function 6(E™ ax — E n ) at some energy E™ ax ~ 0(D m n /T n ) = 
(.D/10 kpc) x 10 9 GeV. Combining these two approximations we obtain 



Normalization to the observed "neutron" excess at ~ 10 9 GeV leads via 
Eq. (|27|) to about 20 antineutrino showers per year |34j . 

A direct v e event in IceCube will make a showering event with poor an- 
gular resolution. Fortunately, neutrino oscillations rescue the signal. Since 
the distance to Cygnus OB2 (D w 1.7 kpc) greatly exceeds the 17 e oscil- 
lation length A sc ~ (-EV/PeV) x 10~ 2 parsecs (taking the solar oscilla- 
tion scale 5m 2 ~ 10~ 5 eV 2 ), the antineutrinos decohere in transit. Replac- 
ing into Eq. ()18|) the most recent SNO result for the solar mixing angle 
9q ~ 32.5° |35| . it is easily seen that the arriving antineutrinos are dis- 
tributed over flavors, with the muon antineutrino flux Fp^ given by the 
factor sin 2 (2 Q )/4 ~ 0.20 times the original Fp e flux. The v T flux is the 
same, and via Eq. (|19j) the v e flux is 0.6 times the original flux. The Cygnus 
region is about 40° below the horizon, hence straightforward calculation 
yields a yearly atmospheric background of 1.5 events. Finally, on the basis 
that HEGRA observations given in Eq. emerge from pion decay, one 
can first estimate the accompanying neutrino flux using Eq. (|25|) . and then 
through Eq. (|27|) verify that the event rate associated with the unidentified 
HEGRA source is even smaller than the atmospheric background. 

In summary, in a few years of observation, IceCube will attain 5a sen- 
sitivity for discovery of the Fe— > n — > V e — > cosmic beam, providing the 
"smoking ice" for the GP neutron hypothesis. 



We have analyzed the possibility of detecting the neutrino counterparts 
of TeV 7-ray observations in a model where 7-rays originate through 7T°- 




(30) 



4. Concluding Remarks 
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decay at the source. We have found that IceCube will attain sensitivity to 
observed neutrinos from the Crab Nebula. The discussion presented here 
can be easily generalized to other sources (see e.g., EH])- Indeed the 
AMANDA-II experiment has already achieved the sensitivity to probe the 
blazar Mrk 501 in its 1997 flaring state if the neutrino and 7-ray fluxes 
are equal jSU- We have also estimated the "essentially guaranteed" V flux 
originated in the decay of neutrons emitted by the Cygnus OB association. 
The expected event rate above 1 TeV at IceCube is dN / dt\ s \ gna \ ~ 20 yr" 1 
antineutrino showers (all flavors) and a 1° directional signal of 4 events, 
well above atmospheric background. 

The detection of (anti)neutrinos pointing towards TeV 7-ray sources or 
to the 10 9 GeV intri guing directional signals would not only provide a final 
answer to the origin of these particles, but also confirm the acceleration of 
protons/nuclei to ultrahigh energies. 
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